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Abstract We report a direct observation of Ge hut nu- 
cleation on Si(OOl) during UHV MBE at 360 °C. Nuclei 
of pyramids and wedges were observed on the wetting 
layer (M x N) patches starting from the coverage of 
> 5.1 A and found to have different structures. Atomic 
models of nuclei of both hut species have been built as 
well as models of the growing clusters. The growth of 
huts of each species has been demonstrated to follow 
generic scenarios. The formation of the second atomic 
layer of a wedge results in rearrangement of its first 
layer. Its ridge structure does not repeat the nucleus. A 
pyramid grows without phase transitions. A structure 
of its vertex copies the nucleus. Transitions between hut 
species turned out to be impossible. The wedges con- 
tain point defects in the upper corners of the triangular 
faces and have preferential growth directions along the 
ridges. The derived structure of the {105} facet follows 
the paired dimer (PD) model. Further growth of hut 
arrays results in domination of wedges, the density of 
pyramids exponentially drops. The second generation of 
huts arises at coverages > 10 A; new huts occupy the 
whole wetting layer at coverages ~ 14 A. Nanocrys- 
talline Ge 2D layer begins forming at coverages > 14 A. 

Keywords Quantum dot arrays ■ Ge hut clusters • 
Nucleation ■ Atomic structure • Molecular beam 
epitaxy • Scanning tunnelling microscopy 
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1 Introduction 

Development of CMOS compatible processes of forma- 
tion of germanium quantum dot (QD) dense arrays on 
the (001) silicon surface as well as multilayer Ge/Si epi- 
taxial heterostructures on their basis is a challenging 
task of great practical significance [l ] [2 l [ l [ l [5 l [6 l [7 l [8 l [9 l 
1101111111211131112] . An important direction of applied re- 
searches in this area is the development of highly ef- 
ficient monolithic far and mid infrared detector arrays 
which could be produced by a standard CMOS tech- 
nology pfTUl[Tll[T2i ri3l[T4]. Such detectors have to com- 
bine high perfection (uniformity, sensitivity, operating 
life, etc.) with high yield and low production price. 
A requirement of CMOS compatibility of technologi- 
cal processes imposes a hard constraint on conditions 
of all phases of the QD array manufacturing starting 
from the stage of preparation of a clean Si surface for 
Ge/Si heterostructure deposition: on the one hand, for- 
mation of a photosensitive layer must be one of the 
latest operations of the whole device production cy- 
cle because otherwise the structure with QDs would 
be destroyed by further high temperature annealings; 
from the other hand, high temperature processes dur- 
ing Ge/Si heterostructure formation on the late phase 
of the detector chip production would certainly wreck 
the readout circuit formed on the crystal. Therefore, 
lowering of the array formation temperature down to 
the values of < 450° C0 is strongly required [TUTT]. and 
the Ge QD arrays meeting this requirement are referred 
to as CMOS compatible ones. 

In addition to the requirement of the low tempera- 
ture of a Ge QD array formation, both high density of 
the germanium nanoclusters (> 10 11 cm~ 2 ) and high 

1 As well as decreasing of the wafer annealing temperatures 
and times during the clean Si(001) surface preparation. 
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uniformity of the cluster shapes and sizes (dispersion 
< 10%) in the arrays are necessary for employment of 
such structures in CMOS IR detectors [12J . The molec- 
ular beam epitaxy (MBE) is known to be the main tech- 
nique of formation of Ge/Si heterostructures with QDs 
[2"1[T5] . A high density of the self-assembled hut clusters 
can be obtained in the MBE process of the Ge/Si(001) 
structure formation when depositing germanium on the 
Si(001) substrate heated to a temperature T gr < 550°C. 
In this case the lower is the temperature of the sili- 
con substrate during the Ge deposition the higher is 
the density of the clusters at the permanent quantity 
of the deposited Ge [TSUI?] . For example, the density 
of the Ge clusters in the array was 6 x 10 11 cm" 2 at 
T gr = 360°C and the effective thickness of the deposited 
germanium layeJl /iq c = 8 A; the cluster density of only 
- 2 x 10 11 cm" 2 was obtained at T gr = 530° C and the 
same value of ho e [IB]- 

There is another approach to obtaining dense clus- 
ter arrays. The authors of Refs. [4"l [H?ll20[|21] reached 
the cluster density of ~ 9 x 10 11 cm" 2 using the pulsed 
irradiation of the substrate by a low-energy Ge + ion 
beam during the MBE growth of the Ge/Si(001) het- 
erostructures at T gr as high as 570°C. 

Obtaining of the arrays of the densely packed Ge 
QDs on the Si(001) surface is an important task but 
the problem of formation of uniform arrays of the Ge 
clusters is much more challenging one. The process of 
Ge/Si(001) heterostructure formation with the Ge QD 
dense arrays and predetermined electrophysical and pho- 
toelectric parameters cannot be developed until both 
of these tasks are solved. The uniformity of the clus- 
ter sizes and shapes in the arrays determines not only 
the widths of the energy spectra of the charge carrier 
bound states in the QD arrays [4] but in a number of 
cases the optical and electrical properties of both the 
arrays themselves and the device structures produced 
on their basis [33) . To find an approach to the improve- 
ment of the Ge QD array uniformity on the Si(001) 
surface it is necessary to carry out a detailed morpho- 
logical investigation of them. 

This article presents the results of our recent inves- 
tigations of several important issues of the Ge dense 
array formation and growth. We have studied the array 
nucleation phase (the transition from 2D growth of the 
wetting layer (WL) to 3D formation of the QD array 
when the nuclei of both species of huts — pyramids and 
wedges [15] — begin to arise on the (M x N) patches 
of WL). We have identified by STM the nuclei of both 
species, determined their atomic structure [TM2"3"] and 

2 I.e. the Ge coverage or, in other words, the thickness of the 
Ge film measured by the graduated in advance film thickness 
monitor with the quartz sensor installed in the MBE chamber. 



observed the moment of appearance the first genera- 
tion of the nuclei. We have investigated with high spa- 
tial resolution the peculiarities of each species of huts 
and their growth and derived their atomic structures 
[231124] . We have concluded that the wedge- like huts 
form due to a phase transition reconstructing the first 
atomic step of the growing cluster when dimer pairs 
of its second atomic layer stack up; the pyramids grow 
without such phase transitions. In addition, we have 
come to conclusion that wedges contain vacancy-type 
defects on the penultimate terraces of their triangular 
facets [23] which may decrease the energy of addition 
of new atoms to these facets and stimulate the quicker 
growth on them than on the trapezoidal ones and rapid 
elongation of wedges. We have shown also comparing 
the structures and growth of pyramids and wedges that 
shape transitions between them are very unlikely [2"31 
124) . Finally, we have explored the array evolution dur- 
ing MBE right up to the end of its life when most of 
clusters coalesce and start forming a nanocrystalline 2D 
layer. 

Below, we present these results in detail. 

2 Methods, equipment and conditions of 
experiments 

The experiments were made using an integrated ultra- 
high vacuum instrument [T8] built on the basis of the 
Riber surface science center with the EVA 32 molecu- 
lar beam epitaxy chamber connected to the STM GPI- 
300 ultrahigh vacuum scanning tunnelling microscope 
[251126 , 27 . This equipment allows us to carry out the 
STM study of samples at any phase of a substrate sur- 
face preparation and MBE growth. The samples can 
be transferred into the STM chamber for the examina- 
tion and moved back into the MBE vessel for further 
processing as many times as required never leaving the 
UHV ambient and preserving the required cleanness for 
STM investigations with atomic resolution and MBE 
growth. 

Initial substrates were 8x8 mm 2 squares cut from 
the specially treated commercial B-doped CZ Si(100) 
wafers (p-type, p — 12 Q cm). After washing and chemi- 
cal treatment following the standard procedure described 
elsewhere |28[|2"5] (which included washing in ethanol, 
etching in the mixture of HNO3 and HF and rinsing in 
the dcionized water), the silicon substrates were mounted 
on the molybdenum STM holders and inflexibly clamped 
with the tantalum fasteners. The STM holders were 
placed in the holders for MBE made of molybdenum 
with tantalum inserts. Then the substrates were loaded 
into the airlock and transferred into the preliminary an- 
nealing chamber where outgassed at the temperature 
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of around 565°C and the pressure of about 5 x 10~ 9 
Torr for about 24 hours. After that the substrates were 
moved for final treatment into the MBE chamber evac- 
uated down to about 10~ n Torr. There were two stages 
of annealing in the process of substrate heating in the 
MBE chamber — at ~ 600°C for ~ 5 minutes and at 
~ 800° C for ~ 3 minutes [TB]. The final annealing at 
the temperature greater than 900°C was carried out for 
nearly 2.5 minutes with the maximum temperature of 
about 925°C (~ 1.5 minutes). Then the temperature 
was rapidly lowered to about 750 °C. The rate of the 
further cooling was around 0.4°C/s that corresponded 
to the "quenching" mode applied in [29]. The pressure 
in the MBE chamber grew to nearly 2 x 10 -9 Torr dur- 
ing the deoxidization process. The surfaces of the silicon 
substrates were completely purified of the oxide film as 
a result of this treatment; more data on the morphol- 
ogy of the prepared Si(001) clean surfaces can be found 
in Refs. [2S1I211I3D] ■ 

Ge was deposited directly on the deoxidized Si(001) 
surface from the source with the electron beam evap- 
oratioi||. The Ge deposition rate was about 0.15 A/s; 
the effective Ge film thickness /iGe was varied from 4 A 
to 15 A for different samples. The deposition rate and 
/iGe were measured by the XTC film thickness monitor 
with the graduated in advance quartz sensor installed 
in the MBE chamber. The substrate temperature T gr 
was 360° C during Ge deposition; the pressure in the 
MBE chamber did not exceed 10~ 9 Torr. The rate of 
the sample cooling down to the room temperature was 
approximately 0.4°C/s after the deposition. 

The samples were heated by Ta radiators from the 
rear side in both preliminary annealing and MBE cham- 
bers. The temperature was monitored with chromel- 
alumel and tungsten-rhenium thermocouples in the pre- 
liminary annealing and MBE chambers, respectively. 
The thermocouples were mounted in vacuum near the 
rear side of the samples and in situ graduated before- 
hand against the IMPAC IS 12-Si pyrometer which mea- 
sured the sample temperature through chamber win- 
dows. The atmosphere composition in the MBE cam- 
ber was monitored using the SRS RGA-200 residual gas 
analyzer before and during the process. 

After Ge deposition and cooling, the prepared sam- 
ples were moved for analysis into the STM chamber in 
which the pressure did not exceed 10 -10 Torr. The STM 
tip was ex situ made of the tungsten wire and cleaned by 
ion bombardment [31] in a special UHV chamber con- 
nected to the STM one. The images were obtained in 
the constant tunneling current (it) mode at the room 
temperature. The STM tip was zero-biased while the 



3 The Si source was switched off during the experiments. 




(a) 



(b) 



Fig. 1 STM image of Ge wetting layer on Si(001): (a) before 
cluster nucleation, h Gc = 4.4 A (U B = -1.86 V, It = 100 pA); 
(b) arising nuclei of pyramidal (1) and wcdgelike (2) huts, /iQe = 
5.1 A (Us = +1.73 V, It = 150 pA). 

sample was positively or negatively biased (U s ) when 
scanned in empty or filled states imaging mode. 

Original firmware [25 , 26 , 27 was used for data ac- 
quisition; the STM images were processed afterwords 
using the WSxM software [32]. 



3 Experimental data and structural models 

3.1 Array and hut cluster nucleation 

Investigating an evolution of the hut arrays we have 
arrived at a conclusion that a moment of an array nu- 
cleation during MBE precedes a moment of formation 
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Fig. 2 STM image of Ge wetting layer on Si(OOl): (a) c(4x2) (c) 
and p(2 X 2) (p) reconstructions within the (M X AT) patches, 
h Gc = 6, A, U s = +1.80 V, It = 80 pA; (b) new formations 
arise on the (M X N) patches due to nucleation of Ge pyramid 
(1) and wedge (2), h Gc = 6, A, U s = +2.60 V, It = 80 pA. 

of the first hut on the WL@ It is not a paradox. Hut 
cluster arrays nucleate when the first hut nuclei arise on 
the (M x N) patch of the wetting layer. This process 
is illustrated by Fig. [TJ An image (a) demonstrates a 
typical STM micrograph of the WL with the (M x N) 
patched structure (/iQe = 4.4 A). This image does not 
demonstrate any feature which might be interpreted 
as a hut nucleus [53]. Such features first arise at the 
coverages ~ 5 A: they are clearly seen in the image 
(b), which demonstrates a moment of the array birth 
(hcc = 5.1 A), and numbered by '1' for the pyramid 
nucleus and '2' for the wedge one (several analogous 
formations can be easily found by the readers on differ- 
ent patches). However, no hut clusters are seen in this 
picture. 

Our interpretation is based on the results reported 
in Ref. [53J which evidenced that there are two different 
types of nuclei on Ge wetting layer which evolve in the 

4 Or, in other words, it foreruns a moment of formation of the 
first {105} faceted cluster with the height-to-width ratio of 1:10 
on WL. 




Fig. 3 Models of nuclei of Ge hut clusters corresponding to the 
images given in Fig. (2{b): (a) a pyramid, (b) a wedge [1 is the 
wetting layer in the plots (a) and (b)]; (c) the models superim- 
posed on the image given in Fig. [2fb), the numbering is the same 
as in Fig. Ob). 



process of Ge deposition to pyramidal and wedge-like 
hut clusters. Having assumed that nuclei emerge on WL 
as combinations of dimer pairs and/or longer chains of 
dimers in epitaxial configuration [33] and correspond 
to the known structure of apexes specific for each hut 
species [T51l2"4"] we have investigated WL patches, 1 mono- 
layer (ML) high formations on them and clusters of dif- 
ferent heights (number of steps) over WL. As a result, 
we succeeded to select two types of formations differ- 
ent in symmetry and satisfying the above requirements, 
which first appear at a coverage of ~ 5 A and then arise 
on WL during the array growth. We have interpreted 
them as hut nuclei, despite their sizes are much less 
than those predicted by the first principle calculations 
[34"] . and traced their evolution to huts. 

The nuclei formation is illustrated by Fig. [5] The 
surface structure of the (M x N) patches is shown in 
the micrograph (a). The letter 'c' indicates the c(4 x 2) 
reconstructed patch, 'p' shows a patch with the p(2 x 2) 
reconstruction [351136] . Both reconstructions are always 
detected simultaneously that means they are very close 
(or degenerate) by energy. The image (b) shows two ad- 
jacent patches reconstructed by the born nuclei: '1' and 
'2' denote the pyramid (a formation resembling a blos- 
som) and wedge nuclei respectively [53]. Their struc- 
tural models derived from many STM images [T51IM1 
|2"3| are presented in Fig. [3Ja, b) and superimposed on 
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Fig. 4 STM image of Ge wetting layer on Si(OOl): (a) h Ge = 
5.4 A {U s = +1.80 V, It = 100 pA) and (b) h Gc = 6.0 A (U s = 
+2.50 V, It = 80 pA). Examples of characteristic features are 
numbered as follows: nuclei of pyramids (1) and wedges (2) [1 
ML high over WL], small pyramids (3) and wedges (2) [2 ML 
high over WL, a T-like wedge 1181 is observed in the image (a)], 
3 ML high pyramids (5) and wedges (6). 



the images of the nuclei in Fig. [2K C )- Note that both 
types of nuclei arise at the same moment of the MBE 
growth. It means that they are degenerate by the for- 
mation energy. An issue why two different structures, 
rather than one, arise to relief the WL strain remains 
open, however. 

It is necessary to remark here that the nuclei are al- 
ways observed to arise on sufficiently large WL patches. 
There must be enough room for a nucleus on a single 
patch. A nucleus cannot be housed on more than one 
patch. So, cluster nucleation is impossible on little (too 
narrow or short) patches (Fig.^b)). 



The hut nucleation goes on during the array fur- 
ther evolution. Fig. H] illustrates this process. An array 
shown in Fig. HJa) (/iGe = 5.4 A) consists of 1-ML nu- 
clei ('1' and '2'), 2-ML and 3-ML pyramids and wedges 
('3' and '5' '4' and '6' mark pyramids and wedges re- 
spectively) Fig. 2Kb) (Iiq c — 6.0 A) demonstrates the 
simultaneous presence of nuclei ('1' and '2') and 2-ML 
huts ('3' and '4') with the growing much higher clusters. 

Hut cluster nucleation on the WL surface continues 
until the final phase of the array life. This peculiarity 
distinguishes low-temperature growth mode from the 
high-temperature one [15] . 



3.2 Structural models 

It is commonly adopted that the hut clusters grow by 
successive filling the (001) terraces of the {105} faces by 
the dimer rows [37]. However, formation of the sets of 
steps and terraces requires the hut base sides to be par- 
allel to the <100> directions. The pyramid nucleus sat- 
isfies this requirement, its sides aline with <I00>. Thus 
the pyramids grow without phase transition when the 
second and subsequent layers are added (Fig. . Only 
nucleus-like structures of their apexes are rotated 90° 
with respect to the rows on previous terraces to form 
the correct epitaxial configuration when the heights are 
increased by 1 ML, but this rotation does not violate the 
symmetry of the previous layers of the cluster. 

A different scenario of growth of the wedge-like clus- 
ters have been observed. Two base sides the wedge nu- 
cleus does not aline with <I00> (Fig. [2Kb)). The ridge 
structure of a wedge is different from the nucleus struc- 
ture presented in Fig. [2Kb) 18,23,21]. It was shown in 
Ref. [23] that the structure of the wedge-like cluster 
arise due to rearrangement of rows of the first layer 
in the process of the second layer formation (Fig. [HKa)). 
The phase transition in the first layer generates the base 
with all sides directed along the <100> axes which is 
necessary to give rise to the {105} faceted cluster. After 
the transition, the elongation of the elementary struc- 
ture is possible only along a single axis which is deter- 
mined by the symmetry (along the arrows in Fig.^a)). 
A formed 2-ML wedge is plotted in Fig. EKb) . A struc- 
ture of the 6-ML wedge appeared as a result of further 
in-height growth is shown in Fig.EKc). The ridge struc- 
tures of the 2-ML and 6-ML wedges is seen to coincide, 
which is not the case for different cluster heights. A 
complete set of the wedge ridges for different cluster 
heights can be obtained by filling the terraces by epi- 
oriented pairs of dimers. 

5 Hereinafter, the cluster heighs are counted from the WL top. 
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Fig. 5 Top views of the pyramidal QDs consisting of (a) 2 and 
(b) 6 monoatomic steps and (001) terraces on the wetting layer 
(1, 2 and 3 designate wetting layer, the first and the second layers 
of the clusters respectively). 



It should be noted also that according to the pro- 
posed model the wedge-like clusters always contain point 
defects on the triangular (short) facets. The defects are 
located in the upper corners of the facets and caused by 
uncertainty of one translation in the position a dimer 
pair which forms the penultimate terrace of the trian- 
gular facet (Fig. [5]). The predicted presence of these 
defects removes the degeneracy of the facets and hence 
an issue of the pyramid symmetry violation which oc- 
curs if the pyramid-to- wedge transition is assumed (this 
issue was discussed in detail in Ref. [IB])- In addition, 
the vacancy- type defects may decrease the energy of ad- 
dition of new atoms to the triangular facets and stimu- 
late the quicker growth on them than on the trapezoidal 
ones and rapid elongation of wedges. These defects are 
absent on the facets of the pyramidal huts. Their tri- 
angular facets are degenerate. Therefore, as it follows 




(a) 




(c) 



Fig. 6 Growth of a wedge-like cluster: (a) reconstruction of 
the first layer of a forming wedge during addition of epi-oriented 
dimer pairs of the second (001) terrace; plots of atomic struc- 
tures of a Ge wedge-shaped hut clusters composed by (b) 2 and 
(c) 6 monoatomic steps and (001) terraces on the wetting layer 
(the numbering is the same as in Fig. [5] d marks a defect arisen 
because of one translation uncertainty of the left dimer pair po- 
sition). 



from our model, the trapezoidal and triangular facets 
of the wedge are not degenerate with respect to one 
another even at very beginning of cluster growth. The 
wedges can easily elongate by growing on the triangular 
facets faster than on trapezoidal ones. Pyramids, hav- 
ing degenerate facets, cannot elongate and grow only in 
height outrunning wedges. This explains greater heights 
of pyramids [TH] ■ 

Analyzing the deduced structural models of pyra- 
mids and wedges, as well as their behaviour during the 
array nucleation and growth, we have come to conclu- 
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sion that shape transitions between the clusters of dif- 
ferent species are prohibited [T51l2"3"ll2"4"] . 

3.3 Facets 

The presented models allowed us to deduce a structure 
of the {105} facets (Fig. EJa)). This model resulting 
from the above simple crystallographic consideration 
corresponds to the paired dimers (PD) [35] rather than 
more recent rebonded step (RS) model [3T)ll4"0"] which is 
now believed to improve the previous PD model by Mo 
et al. 

A direct STM exploration of the {105} facets con- 
firms the derived model. Being superposed with the 
empty state STM image of the cluster {105} facet it 
demonstrates an excellent agreement with the experi- 
ment (Fig.EJb)). A typical STM image of the QD facet 
is presented in Fig. [5J Characteristic distances on the 
facets are as follows: ~ 10.5 A in the <100> directions 
(along the corresponding side of the base) and ~ 14 A 
in the normal (<051>) directions. The facets are com- 
posed by structural units which are outlined by ellipses 
in Fig. [5Ja) and can be arranged along either [110] or 
[110] direction on the (001) plane. We have interpreted 
them as pairs of dimers. Their positional relationship 
is obviously seen in the 3D micrograph presented in 
Fig. lib). 

Dangling bonds of the derived {105}-PD facets, due 
to high chemical activity, may stimulate Ge atom ad- 
dition and cluster growth. Thus less stability and higher 
activity of the {105}-PD facets compared to the 
Ge(105)/Si(105)-RS plane, which is usually adopted in 
the literature for simulation of hut {105} facets, may 
cause fast completion of hut terraces during epitaxy 
and be responsible (or even be necessary) for hut for- 
mation and growth. 

3.4 Cluster density and fractions 

Fig- Et a ) plots the dependence of the cluster density on 
ft-Gc for different clusters in the arrays. It is seen that 
the density of wedges rises starting from D w ps 1,8 x 
10 11 cm~ 2 at the beginning of the three-dimensional 
growth of Ge (the estimate is obtained by data ex- 
trapolation to 1iq c = 5 A) and reaches the maximum 
of ~ 5 x 10 11 cm~ 2 at 1iq c ~ 8 A, the total den- 
sity of clusters at this point ~ 6 x 10 11 cm~ 2 
is also maximum. Then both D w and slowly go 
down until the two-dimensional growth of Ge starts at 
h Gc ~ 14 A and D s w D w - 2 x 10 11 cm -2 (the con- 
tribution of pyramids D p to Dy, becomes negligible — 
~3x 10 10 cm -2 — at this value of Hq c ). The pyramid 



'<I00> 




(a) 



Fig. 7 (a) A structural model of the {105} facet of hut clus- 
ters derived from the plots given in Figs. [5] and [6] corresponds 
to the PD (pairs of dimers) model I38| . Sa and Sb are com- 
monly adopted designations of the monoatomic steps 41 , atoms 
situated on higher terraces are shown by larger circles, (b) 
The schematic of the facet superimposed on its STM image 
(4.3 X 4.4 nm, U s = +3.0 V, It = 100 pA), the [100] direction 
is parallel to the corresponding base side, the steps rise from the 
lower right to the upper left corner. 

density exponentially drops as the value of ft-Gc grows 
(Dp s» 5 x 10 u cxp{-2.0 x 10 7 h Gc }, h Gc is measured 
in centimeters). The maximum value of D p ks 1.8 X 
10 11 cm -2 obtained from extrapolation to Iiq c = 5 A 
coincides with the estimated initial value of D w . 

The graphs of cluster fractions in the arrays ver- 
sus ha,, are presented in Fig. EJb). Portions of pyra- 
mids and wedges initially very close (~ 50% at ft,Qe ~ 
5 A) rapidly become different as Hq c rises. The con- 
tent of pyramids monotonically falls. The fraction of the 
wedge-like clusters is approximately 57% at the early 
stage of the array growth (/iQ e = 6 A) and becomes 
82% at hoe = 8 A. At further growth of the array, the 
content of the wedges reaches the saturation at the level 
of approximately 88% at hc c = 10 A. 

The inference may be made from this observation 
that contrary to the intuitively expected from the con- 
sideration of symmetry, the wedge-like shape of the 
clusters is energetically more advantageous than the 
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Fig. 8 (a) 2D and (b) 3D STM images of the same area on 
Ge hut cluster facet (h Go = 10 A, T gr = 360°C, U B = +2.1 V, 
It = 80 pA). The sides of the cluster base lie along the [100] 
direction; structural units revealed on the free surfaces of the 
(001) terraces and interpreted as paired dimers are marked out. 



| 

d. 10 : 



5 6 7 8 9 10 11 
Ace (A) 



(a) 



5 6 7 8 9 10 11 12 13 14 15 
A<;,(A) 



(b) 



Fig. 9 (a) Density and (b) fraction of the Ge clusters in the ar- 
rays formed at T gr = 360°C (□ marks the pyramids, ■ designates 
the wedges, O is the total density). 



pyramidal one, and the more advantageous the more Ge 
atoms (and the more the number of terraces) constitute 
the cluster. The probability of nucleation appears to be 
close to 1/2 for both wedges- like and pyramidal clus- 
ters at the initial stage of the array formation and low 
growth temperatures. Then, as the array grows, the for- 
mation of pyramids becomes hardly probable and most 
of them, which have already formed, vanish whereas 
the nucleation and further growth of wedges continues. 
The Ge pyramids on the Si(001) surface turned out to 
be less stable objects than the wedges. 

Notice also that at T gr = 360° C and the flux of Ge 
atoms d/iQe/di = 0.15 A/s, the point /iqo = 10 A is 
particular. Not only the fraction of pyramids saturates 
at this point but the array in whole has the most uni- 
form sizes of the clusters composing it (Tigs. HI [TUland 
[TTj) . This is concluded by us not only on the basis of 
analysis of the STM images of the Ge/Si(001) arrays 
but also from the data of the Raman scattering by the 
Ge/Si heterostructures with different low-temperature 
arrays of Ge quantum dots H2"ll43] . We refer to such 
arrays as optimal. 




X: 51.6nm 



X: 99.7nm 



Fig. 10 STM 2D and 3D micrographs of Ge hut cluster dense 
arrays at different coverages (T gr = 360° C): (a), (b) h Gc = 8 A 
[(a) 50.6 X 49.9 nm, w is the wetting layer, U s = +2.0 V, It = 
80 pA; (b) U s = +2.0 V, It = 100 pA]; (c),(d) h Gc = 10 A 
[(c), (d) U s = +2.1 V, h = 100 pA]. 




X: 32.9nm 




X: 55.6nm 
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Fig. 11 STM topographs of Ge hut cluster dense arrays at 
different coverages (T gr = 360°C): (a), (b) h Gc = 14 A [(a) U s = 
+1.75 V, h = 80pA, (b) U s = +3.0 V, It = 100pA]; (c), (d) 
h Ga = 15 A [(c) Us = +2.0 V, h = 120 pA, (d) 20.3 x 20.4 nm, 
Us = +3.6 V, It = 120 pA], w indicates the wetting layer patches, 
i shows a distorted small Gc island 3 ML high over WL. 
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3.5 Array life cycle 

A qualitative model accounting for the presence of the 
particular point at the low-temperature array growth is 
simple. The case is that at low enough temperatures of 
the array growth, the new Ge cluster nucleation com- 
petes with the process of growth of earlier formed clus- 
ters. The height of the dominating wedge- like clusters 
is observed to be limited by some value depending on 
2gr@ At small hoe, Ge clusters are small enough and 
the distances between them are large enough compared 
to the Ge atom (or dimer) diffusion (migration) length 
on the surface for nucleation of new clusters on the Ge 
wetting layer in the space between the clusters (Figs.|H 
fTOfal.fb)). At h Gc = 10 A and the above dh Gc /dt val- 
ues, the equilibrium of parameters (cluster sizes and 
distances between them, diffusion length at given tem- 
perature, Ge deposition rate, etc.) sets in, the rate of 
new cluster nucleation is decreased and the abundant 
Ge atoms are mainly spent to the growth of the avail- 
able clusters (Fig. [EIc),(d)). After the clusters reach 
their height limit and in spite of it, Ge atoms continue 
to form up their facets. As soon as most of the clus- 
ters reach the height limit, nucleation of new clusters 
becomes energetically advantageous again and the nu- 
cleation rate rises. The second phase of clusters appears 
on the wetting layer and fills whole its free surface as 
fiQe is increased (Fig. [11]). Further increase of hc c re- 
sults in two-dimensional growth mode. It is clear now 
why the array is the most homogeneous (optimal) at 
T gr = 360°C and }iq c — 10 A whereas the dispersion of 
the cluster sizes is increased at higher and lower values 
of hcc because of the small clusters containing in the 
array. It is clear also that the optimal array will ap- 
pear at different value of h Gc when T gr or dhce/dt are 
different. 

As it follows from the data presented in this section 
and Section l3~Tl the Ge hut array evolution and life cycle 
goes through three main phases: at T gr = 360°C, the 
array nucleates at h Gc ~ 5 A (Fig.[T]), it reaches ripeness 
and optimum to 1iq c ~ 10 A (Fig. [TU)) and finishes its 
evolution at Hq c ~ 14 A by filling whole the surface 
(Fig. IT!"]). Most of clusters start coalescing (Fig. Qjjb)) 
and 2D growth begins at greater ho c (Fig. [TTJc)). 

Nevertheless, free areas of WL still remain even at 
fiQ C — 15 A (Fig. ITj7d) ). The structure of the parches 
('u>') stays the same as in the beginning of the array 
formation although the WL regions are surrounded by 
large huts. Small 3D islands ('i'), although very dis- 



6 Note that we did not observe a height limitation of pyramids. 
We suppose that they may give rise to one of the types of array 
defects — huge clusters — which sometimes appear among ordinary 

huts [3D. 



torted, are still recognizable on WL between the large 
huts. The hut nucleation on WL goes on even at as 
high coverages as 15 A when virtually total coalescence 
of the mature huts have already happened. 

4 Conclusion 

In summary, we have studied the array nucleation phase 
and identified the nuclei of both hut species, determined 
their atomic structure and observed the moment of ap- 
pearance of the first generation of the nuclei on WL. 
We have investigated with high spatial resolution the 
peculiarities of each species of huts and their growth 
and derived their atomic structures. We have concluded 
that the wedge- like huts form due to a phase transition 
reconstructing the first atomic step of the growing clus- 
ter when dimer pairs of its second atomic layer stack 
up; the pyramids grow without phase transitions. In 
addition, we have come to conclusion that wedges con- 
tain vacancy-type defects on the penultimate terraces 
of their triangular facets which may decrease the energy 
of addition of new atoms to these facets and stimulate 
the quicker growth on them than on the trapezoidal 
ones and rapid elongation of wedges. We have shown 
also comparing the structures and growth of pyramids 
and wedges that shape transitions between them are 
impossible. And finally, we have explored the array evo- 
lution during MBE right up to the concluding phase of 
its life when most clusters coalesce and start forming a 
nanocrystalline 2D layer. 
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